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Abstract 
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We consider a spinor Bose-Einstein condensate in its polar ground state. We 
analyze magnetization waves of a finite amplitude and show that their nonlinear cou- 
pling to the density waves change the dependence of the frequency on the wavenum- 
ber dramatically. In contrary, the density wave propagation is much less modified 
by the nonlinearity effects. A similar phenomenon in a miscible two-component 
condensate is studied, too. 

PACS number: 03.75.Fi 

Recent advances in experimental creation of multicomponent atomic Bose-Einstein 
condensates |1|, |2|, |3| have given rise to an interest to physical properties of such systems. 
There are numerous works on the properties of degenerate Bose gas mixtures in magnetic 
traps related to both the ground state || and the collective excitations ||. In Ref. |5[ 
the early work || related to a homogeneous Bose gas mixture is generalized to the case 
of presence of external harmonic trap potential. Of course, the number of branches 
of the dispersion law is equal to the number of different components in a mixed BEC 
Due to non-zero interaction between them, normal mode oscillations imply simultaneous 
mutually coherent motion of the components. In the present paper, we, however, consider 
firstly a multicomponent BEC of another kind, namely, a spinor BEC. Such a degenerate 
j> \ quantum system can be created in an optical trap, where all the atoms are confined 

practically independently on rrif, their momentum projection to any arbitrary axis. Such 
an independence of confinement on the spin orientation is a striking feature and a key 
advantage of an optical trap, well justified experimentally 0, |3|. Under such a condition, 
the spin orientation becomes a new degree of freedom. The differences and similarities 
between a two-component BEC with fixed m/'s for both the components and a spinor 
BEC in context of our study will be discussed later, at the end of this paper. 

Now we have to note that in all the cited works on collective excitations in multi- 
component BECs as well as in the seminal works on spinor BEC dynamics oscillation 
amplitudes were assumed to be small enough to provide linearization of the set of coupled 
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time-dependent Gross-Pitaevskii equations (GPEs). Then a proper linear transformation 
yields the equations of the harmonic oscillator type for the normal modes. However, the 
GPE is essentially nonlinear, thus the effects of finite amplitude of oscillations take place. 
There are some approaches to take nonlinearity into account. The first one is to find 
particular solutions of the GPE in a form of solitons (see, e.g., the recent work |J and ref- 
erences therein). The second one is to find oscillating nonlinear solutions those in the case 
of infinitesimally small oscillation amplitude coincide with corresponding eigenfunctions 
of the linearized version of the GPE or of the equivalent set of quantum hydrodynamical 
equations. An elegant formalism has been developed for nonlinear oscillations of a scalar 
BEC in a harmonic trap in the Thomas-Fermi regime |J. It has been found that the 
nonlinear effects become important, if the fraction of mass of a scalar BEC involved into 
oscillatory motion is comparable to unity. 

In the present paper, we study validity of the approximation based on linearization of 
the GPE by proceeding in the following way. We consider plane waves in a spatially ho- 
mogeneous multicomponent BEC. This can serve as a WKB approximation for excitations 
in a trapped BEC, if the excitation wavelength is much smaller than the atomic cloud 
size. Moreover, such an approach allows us to use in the most direct and straightforward 
way the standard technique of expanding a solution in series in a certain small parameter, 



known as a standard perturbation theory in classical mechanics [10]. Also consideration 
of plane waves in a translationally invariant BEC provides a possibility of comparing the 
results to the strict analytic formulae of Refs.[|. [/]]. 

The main result of our work is that certain modes in a multicomponent BEC exhibit 
strongly nonlinear behaviour, namely, even for a relatively small wave amplitude the 
effects of anharmonicity become significant. This effect is absent in a scalar BEC case. 

Let us consider a spinor BEC composed of atoms with the spin / = 1 at zero temper- 
ature. The GPE governing evolution of the complex order parameter (macroscopic wave 
function) ?/>(r, t) of the BEC reads in the mean field approximation as follows 0: 

in ^fi = 'Im 7 *^ ~^ + hc °^*^ + MW) ■ (*V0, (1) 

where f is the single-atom angular momentum operator, a vector with Cartesian com- 
ponents being 3x3 matrices, M is the mass of an atom, \x is the chemical potential. 
Interaction constants defined as hc = (go + 2g 2 )/3, hc 2 = (g 2 — <7o)/3, Qf = 4:Trh 2 a F /M, 
and ap is the s-wave scattering length for a pair of slow atoms with the total angular 
momentum F equal to or 2, respectively. Practically, the magnitudes of these two scat- 
tering lengths are close each to other, so |c2/cq| <C 1. The order parameter if> has three 
components, corresponding to the momentum projection to the 2-axis rrif = 0, ±1: 
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where n is the total equilibrium density of the BEC. if) means the transposed vector. In 
other words, the ground state components of the vector ( are normalized by the condition 

E |d 9 i oW)2 = i. (2) 



-1 



We assume the repulsive interaction of atoms in the BEC, i.e., cq > 0. For the sake 
of definiteness, we assume also that c 2 > 0. It follows from the latter condition that the 
ground state of such a system is a so-called polar state J7J. This means that, in the mean 
field picture, all the atoms have zero momentum projection to a certain axis. This state 
is degenerate with respect to orientation of this axis. Let this axis be the z-axis, so in the 
equilibrium, when the time derivative of ip in Eq.([l]) is equal to zero, 

>.(ground) ^ ^(ground) -. 

S>±1 — U > SO = *■■ 

The chemical potential of the BEC in the polar state is \i — Cr,n. 

Before writing Eq. ([!]) in explicit form, we introduce the new unknown functions: £± = 
(Ci ± C*i)/v2, r] p = Re Co — 1, Vi = Im Co- Then Eq.(|I|) can be transformed to the set of 
equations 

<W(£+C - ££- + 2th + 2w)C- + 2(2^ p + $£+], (3) 

c 2 n[(e+c - e;e- + ^ + 2 W )£+ + Ht-i w 

d h 

-Q^i = -^^ 2 Vp + ^ nri P + Con[(C + ^+ + C^- + ^Vp + V 2 p + Vi)Vp + 

^e + + ££- + V 2 } + c 2 n[2C + Z + V P + 2££ + + m- + Z+Ovl (5) 

c 2 n[2e^-^ + (e;e_ + e+e)^ + ae- + e+ei- (6) 

If we neglect all the nonlinear terms in Eqs.(|3] — §), then we get immediately the solutions 
in the form of plane monochromatic waves and the corresponding dispersion laws ]7| . The 
first mode is the density wave, it corresponds, in the linear approximation, to the periodic 
oscillations of the m/ = component of the order parameter only (i.e., of rj p , rji), while 
£ + and £_ remain zero. Density waves in a spinor BEC are the same as sound waves in a 
scalar BEC. The dependence of the frequency u^o of the density waves on the wavenumber 
k is of the Bogoliubov's type, uJ%(k) = uj r {k)[uj r {k) + 2c n], where ui r {k) = hk 2 /(2M) is 
the recoil frequency associated with the kinetic momentum Hk. Another branch of the 
excitation spectrum in a spinor BEC is related to magnetization waves. Left and right 
circularly polarized magnetization modes are degenerate and in the linear regime their 
frequency is given by the formula u^ nQ (k) = u r (k)[u r (k) + 2c2n]. The quantum mechanical 
mean values of the atomic magnetic momentum operator are proportional to the £ + and 
£+ for the left and right polarization, respectively. 

Now we can determine the effects of nonlinearity on magnetization wave propagation, 
using the perturbation theory of classical mechanics Ql0j| . Namely, we expand our unknown 



functions in series: £ + = X^o£+ ' where £+ is proportional to the j-th power of a 
certain small parameter e (in fact, the square of the magnetization amplitude can be 
naturally regarded as this parameter). Similar expansions hold for the remaining three 
functions. The zeroth order approximation can be taken also in the form of the plane wave, 
£+ = A + sin (tot— kr), but with the frequency to shifted with respect to the non-perturbed 
value uJ m Q. The validity of this method is restricted to the case of small resulting correction 
to the frequency, \(uj — u m0 )/uj m0 \ <C 1. Also we take £_ = u> r (k)~ 1 u>A + cos (cut — kr), 
r)^ = 0, 7]\ = 0. The difference between uj and uj m0 can be also represented as a series 
in e, beginning from the term of order of e 1 . 

To find the correction to the frequency of a magnetization wave, we make the following 
transformation of our set of GPEs. We add to and substract from the right-hand side of 
Eq.(|3p the term uj 2 ^ + /oj r (k). Then we note that our zeroth order approximation satisfies 
the set of equations —d^/dt = uj 2 t; + /uj r (k), d£ + /dt = uj r (k)^_ identically. The remaining 
terms must be regarded as a perturbation leading to the frequency shift in higher orders 
of approximation. Eqs.([3] — 16]) must be satisfied in every order in e separately, i.e., one 
must group all the terms of order of £ J in the right-hand side and equalize them to the 
0(si) part of the left-hand side of the equation. We restrict our analysis to the linear 
order in e when we obtain 
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JX = uj r (k)C-' H tttCo^ sin (ut-kr)+ 9/ ,.. cos (wt-kr 

at u r (k) L ^r(^) 

A^cos(cjt-kr). (8) 

Here the symbol { ... }^ means that only linear in e ~ A 2 + contribution to the expression 
in the curly brackets is retained. A + is taken to be real, without loss of generality. 
Eqs.(0, H) can be easily reduced to the following differential equation 
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A+ sin (out - k r) + CA\ sin [3(ut - k r)] = 0. (9) 

Here C is a certain combination of various frequency parameters of the problem; its 
calculation is not needed for determination of the correction to the wave frequency in the 
lowest order. 

Eq.(|9|) is inhomogeneous, and the presence of resonant source term proportional to 
sin (cut — k r) leads to occurrence of oscillations with amplitude growing linearly in time in 
the solution for £+ . And the very essence of the method used here |1(J is avoiding of these 
non-physical (secular) solutions proportional to tsm(oot — kr) by setting the prefactor of 
the resonant term to zero. Thus, to the lowest order in square of the wave amplitude, the 
magnetization wave frequency is given by the expression 
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In the two limiting cases (of the short and long wavelength) we obtain 

u 2 = w 2(jfe) + -u 2 k 2 A 2 + , hk > Mu 2 , (11) 



uji = ujk 2 + 6ulk 2 (^Y A 2 + , hk<^Mu 2 , (12) 
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respectively. Here up = JTicfti/M are the velocities of propagation of density (F = 0) 
and magnetization (F = 2) waves of infinitely small amplitude in the long wavelength 
limit. So we can conclude that the effects of nonlinearity are small until 

hk > Mu , (13) 

Mu 2 <hk< Mu , (14) 




hk < Mu 2 . (15) 

This is interesting that the condition of the nonlinearity smallness coincides with the 
trivial condition of smallness of A 2 + in comparison to the sum of squares of absolute 
values of all the three Cm™ 1 " 1 ^ i n the ground state which is unity, according to Eq.(0), 



only in the short wavelength limit of Eq. fllTf) . In the other other cases [Eqs.(|l4], 15)], even 
small but finite excitation amplitude can result in a significant modification of the wave 
propagation. 

It is easy to show that Eqs.(|5|, §), in the case of magnetization waves, have no resonant 
terms in their right-hand sides in the first order in e and hence do not contribute to the 
evaluation of the corresponding correction to the wave frequency. 

Density waves can be analyzed in the similar way, and the lowest-order correction 

results in the formula 

3 

^ 2 = ^do + -^c nuj r (k)A 2 p , (16) 

where A p is the amplitude of oscillations of r] p . For all the momenta k, the correction is 
small provided that A p <C 1, i.e., nonlinear effects play less role for waves of this type, 
in contrary to magnetization waves. Eq.flTo]) also applies to sound waves in a single- 
component (scalar) BEC 

Since u> m0 does not depend on cq but the latter quantity appears in the right-hand 
side of Eq. (|i~0|) , we conclude that nonlinear coupling to density waves plays the key role in 
modification of the magnetization wave frequency. In opposite, Eq.flTB]) does not contain 
c 2 , so a traveling density wave is not coupled to magnetization modes. 

Now let us discuss briefly the case of a mixture of two BECs each having the fixed value 
of rrif or, equivalently, of two scalar BECs. Here we need first to introduce the coupling 
constants g^j = 2nhajij(Mj + Mj>)/(MjMj>), where Mj is the mass of an atom of the 
j-th kind, cijij is the s-wave scattering length for a pair of atoms of j-th and j'-th kind, 
j', j = 1, 2. The dispersion laws for the two excitation branches were obtained in the 
analytic form in Ref. (see also Ref.||). If all the three relevant scattering lengths are 
positive, the criterion of stability of a homogeneous BEC mixture against phase separation 



is simply g\ 2 < JQwQii- In this case, the eigenmode frequencies are positive for all the 
values of the momentum k. For the sake of simplicity, we consider in our paper the 
case of equal atomic masses, Mi = M 2 = M. Then the eigenfrequencies are simply 



J± = u r (k)[u r (k) + 2A ± ], where A± = [g n ni + g 22 n 2 ± \J(gnni - g 22 n 2 ) 2 + 4#? 2 nin 2 ]/2, 
Hi, n 2 are the equilibrium number densities of the components, u r (k) is the same quantity 
as defined above. 

The order parameter perturbation for the j-th component reads as 5ipj = ^fn{A,j [sin(u;£- 
kr) + iu^ 1 {k)oj cosiyjt — kr)]. After some tedious but straightforward calculations, anal- 
ogous to those described above and valid under the same condition of smallness of the 
frequency correction, we arrive at the following formula for the shifted, due to the non- 
linearity effects, wave frequency: 
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Here the upper sign corresponds to the case of B + 7^ 0, 5_ = 0, and the lower sign 
corresponds to the opposite case, B + = 0, £>_ 7^ 0. Here the eigenmode amplitudes are 
defined as 

B + = cos6 g A 1 + J— sin 6 g A 2 , B„ = -sind g Ai + J— cos6 g A 2 . (18) 

Also we set by definition 

g + = gu cos 4 8 g + 2g 12 cos 2 9 g sin 2 6 g + g 22 sin 4 6 g , (19) 

g- = gu sin 4 9 g + g 22 cos 4 9 g , (20) 



g 22 n 2 - guru + J(g 22 n 2 - gun^ 2 + m 2 n x n 2 

tan^ = s - — = . (21) 
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Eq. (|17D is of the form similar to Eq. (|10D and leads to similar restriction on the wave 
amplitude. If we the two BECs are composed of atoms accumulated on two different 
magnetic or hyperfine sublevels of the ground internal state, the difference between g\ 2 
an d yjQwQii is relatively small, and the lower-frequency mode is extremely sensitive to 
the effects of nonlinearity in the long wavelength limit. One should note that both the 
branches of excitation spectrum of a two-component BEC in an external magnetic field 
are sensitive to nonlinear effects for small fc's, while the spinor BEC collective excitations 
exhibit different behaviour: the nonlinearity effects are important for magnetization waves 
much more than for density waves. 

In summary, we should note that the studied effects of nonlinearity in wave propa- 
gation in a BEC are related to the Beliaev damping [JTTJ] (cf. the closely related recent 
publication |12] on an efficient damping of a relative motion of two condensates in a trap 
by nonlinear interaction). The Beliaev damping is also described the cubic nonlinear term 
in the GPE. It is, in fact, decay of a quantum of collective excitation to two quata of lower 
energies, provided that the energy and momentum are conserved. This process results 
in occurrence of an imaginary part of the wave frequency (the damping constant). We 
calculate in the present paper the real small addend to the wave frequency. While the 



Beliaev damping becomes less important when k approaches zero, the nonlinear correc- 
tions to the magnetization mode in a spinor BEC and to both of the modes in an usual 
two-component BEC become more pronounced. 

Finally, we present a numerical example. The ground state of a spinor BEC of sodium 
atoms with / = 1 is just a polar (antiferromagnetic) state 0. We take (a +2a 2 )/3 ~ 5 nm, 
(02 — a )/3 ~ 0.08 nm and set n rs 10 14 cm -3 . Let the excitation wavenumber be of about 
3.5 ■ 10 3 cm -1 (the corresponding wavelength is several times smaller than the atomic 
cloud size in an experiment with large number of atoms in a trap like that of Ref . , so 
the WKB approximation still remains satisfactory). For A + — > the linear theory |/J 
gives the magnetization wave frequency u> m o « 300 s _1 . However, if A + « 0.044, in other 
words, only [1 + uj~ 2 [k)uj 2 ]A 2 + / 2 « 0.005 of the total mass of the BEC is involved into the 
motion, then the frequency rises by one-third of its primary value and becomes equal to 
400 s _1 , according to Eq. (|T0"D . Similarly, one can expect a strongly nonlinear behaviour of 
low-lying magnetization modes of a spinor BEC in an optical trap of a finite size, since the 
trapped BEC spectrum ought to reveal the most important qualitative features present 
in the translationally invariant case, as it has been shown for two-component BECs in 
magnetic traps ||. 
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